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1 Introduction 



The possibility that the fundamental constants of nature might vary with time has been an 
object of speculations for many years In a modern language, however, one should say that 
the prospect of a time variation in the vacuum expectation value (hereafter vev) of the Higgs field 
seems more plausible than the time variation of the Fermi coupling constant Gp or of the electron 
mass me 0. 

Changing the vev of the Higgs field has many physical effects, four of which have astrophysical 
consequences: Gp changes, the electron mass rrie changes and the nuclear masses and binding 
energies change. All of these effects alter Big Bang nucleosynthesis. On the other end the change 
in rrie is the only effect relevant for the cosmic microwave background spectrum 0. 

How much do we know about possible variations of rrie? We should remember that one of 
the recurring themes in the physics behind the fundamental constants is that their values are 
rarely determined by a direct measurement. The example of the electron mass illustrates how 
the information that leads to the values of the constants can be indirect and how different paths 
provide redundant constraints on their values. To obtain the best values, all of this information 
are taken into account simultaneously; in the approach of the 1998 adjustment 0], the information 
is divided into three categories: input data, observational equations, and adjusted constants. The 
observational equations are theoretical expressions that give values of the quantities in the input 
data category as functions of the adjusted constants. The adjusted constants are a suitably chosen 
set of fundamental constants that are determined by the adjustment. The adjustment's role is to 
find the values that best reproduce the input data by means of the theoretical expressions. 

There are recent studies that show, for instance, how a change in me alters the CMB fluctuation 
spectrum. There it is assumed that the variation in me is sufficiently small during the process of 
recombination so that, one needs only consider the difference between me at recombination and 
me today, see and also 0. Furthermore, it has been pointed out that MAP and PLANK 
experiments might be sensitive to variations as small as |Ame/me| ~ 10^^ — 10^^. 

Although all of these considerations look very promising, we are still missing an important in- 
gredient in the discussion: can a time variation in the electron mass (any mass) be made formally 
consistent with the Standard Model of strong and electroweak interactions? In other words, do 
we have to assume an ad hoc time dependence in these parameters or do we have some explicit 
time variation which is consistent with the mathematical structure of the Standard Model? Fur- 
thermore, a time-dependent mass is a ill-defined concept since there is no stationary state in a 
time-dependent external field. 

In this paper we show that the classical equation of motion for the Higgs field in the Minimal 
Standard Model admits a time- dependent solution which can be given in terms of a Jacobian el- 
liptic function, therefore a doubly-periodic function of time. This solution is explicitly constructed 
in Sect. 2. In Sect. 3 we consider the coupling of fermionic fields with the Higgs sector and study 
the effect of their propagation in a time-dependent Higgs background. The three-momentum p 
is conserved since the background depends only on time and we can factorize the usual term, 
exp(i p • x). In the limit of very large primitive period of the Higgs background the Dirac equation 
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for the electron field can be solved in WKBJ approximation. The solution can be cast into the 
form of a plane-wave with a time-dependent distortion factor which becomes arbitrarily small 
exactly in the limit of infinite primitive period of the Higgs background. 



2 A time dependent solution for the Higgs vev 

The relevant part of the Standard Model Lagrangian that we are interested in is 

= -d.K^d^K -^^K^K-^\ {k^k)\ (1) 
where is a complex iso-doublet, 

^=^(^ + #"0(j), (2) 
and where we have neglected gauge couplings. We introduce two new quantities as follows: 

Ml=Aj^M\ ^^ = p-^Ml, (3) 

where g is the SU (2) coupling constant and, for a constant vev, is the bare Higgs boson mass. 
As usual we perform a shift in the ■j/'-field 

M 

^ = 2 — ^ + H. (4) 
9 

The parameter (3 will be adjusted in perturbation theory so that, order-by-order, < 0|if|0 >= 0. 
The well known solution which is at the basis of the so-called spontaneous symmetry breaking 
mechanism is $^ = 1. For a constant, non-zero, value of < Q\ip\Q > M is the bare VT-boson mass 
and the bare Higgs boson mass. In the present case we allow $ to be function of time and 
arrive at the following equation, where according to the usual procedure /? = at tree level: 

^$_iM2$ fl-<l>2) =0. (5) 
A solution to this equation is given in terms of Jacobian elliptic functions 0, 

N', = % Ll = 2e-1, H, = ^, (6) 
and k is any positive real number. Note that we adopt the definition 

f4> f \ —1/2 

u = y dt [X — k^^ sin^tj , = am(u), 

cn(M) = cn(M, A;) = cos(am(u)), sn(u) = sn(u. A;) = sin(am(u)), 

r 1 1/2 

dn(M) = dn(u, fc) = 1 — A;^ sin^(am(u)) . (7) 
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Such solutions are well known in the literature for both (l + l)-dimensional and (3 + l)-dimensional 
variants of 0^ theory ||^. The functions cn{z, k), dn{z, k) have the following properties: 



periods : 4 K,2K + 2iK' 2K,4iK' 

zeros: (2m + l)K + 2raK' (2 m + 1) K + (2 n + 1) i K' (8) 

poles : I3^n = 2 mK + (2 n + l)i K' 

where n,m are integers and K(A;) is the complete elliptic integral of the first kind 

K{k) = F{^,k), K'{k) = F{^,k% (9) 

with /c' = (1 — /c^)^/^. Clearly we are interested in solutions with large values of k since in that 
limit 

cn{Hk,t,k) = dn(^T,^^ r^l-^sin^T + ol^^y T = kHkt, (10) 

and the vev of the Higgs field is approximately constant with time periodic fluctuations suppressed 
by a factor k^"^. The Lagrangian becomes 

= ~ d,Hd,H + 2 ^ [$2 - (/3 - i M^) $2 - i Ml + O [h^] . (11) 
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This result, and the vev of the if-field in higher orders of perturbation theory require some 
additional comment. There is a total derivative to be considered. 

Therefore, to quantize the ff-field, we need to consider a time interval {— t^, +1"^} 

r^ = (2m + 3)K, r«=(2m+l)K, (13) 

where m is an integer, 4 K is the primitive period of $ and $(— r^) = $(Tfl) = 0. For the quantum 
fluctuations H we impose periodic boundary conditions. At the end the limit m ^ oo will be 
taken. 

In this paper, however, we are not so much interested in the Higgs sector, i.e. Higgs mass and 
Higgs self-couplings, or in the problem of time dependence of the cosmological constant but rather 
we concentrate on the effect of a time-depended Higgs vev on the masses of elementary particles. 
For this reason we will analyze in the following section the Higgs-fermion Yukawa couplings. 

3 Higgs-fermion interaction 

The relevant piece of the Lagrangian for arbitrary d fermion fields will be 



where ipL ^ left-handed doublet and u{d)^ are right-handed singlets of SU{2) with being the 
charge-conjugate of K. For $ constant are just the bare up, down masses. For the (z/g, e) 

doublet we find ^ 

Cm = -T^epJ^e-epe - -gj^ijj*ee. (15) 

Here me is an arbitrary free parameter with dimension of mass which cannot be identified with 
the bare electron mass since the fermion fields are moving in a time-dependent background. The 
neutrino remains decoupled from the Higgs field while the Dirac equation for the electron becomes 

(^ + me$)e = 0. (16) 

If we split the field e into upper/lower two-dimensional components, e+/e_, we obtain the following 
equations: 

iXa''daex + i\dte_x + m,^e^x = 0, (17) 
where a", a = 1, 2, 3 are Pauli matrices and A = ±1. To find a solution we introduce 

eA = e^^"^''X., (18) 
and derive the corresponding equation for x^, 

(T''paX,-idtX-,- ^me<l>X-,=0. (19) 
The above equation is solved by introducing a new set of variables, scalar and vector modes: 

X, = ^iSx + ^Vx^aa) (j) . (20) 

Next we introduce a set of matrices (J°'{p), such that 

P± = ^ (PxTiPy) , P = IpI, 

Uu = N, 1 = 1,2, Uui2i) = ±V2N^^, N-^ = 2-^. (21) 

P + Pz P + Pz 



With this definition it is easy to prove that the matrices 

a%p) = U^a^U, (22) 

satisfy the following properties: 

a3(p) = a-p, ai'2(p) = a-ei,2, (23) 

with 63 = p/p and 



4 



Using these matrices we derive a useful decomposition for the vector modes, 

i=l,2 



1=1,2 

By equating the coefficients of 1 and cr(p) we obtain two separate systems of equations 
the SL (scalar-longitudinal) and T (transverse) modes: 

^ ^(A) + $ S{\) - p K(-A) = 0, 
at 



dt 
d 

di 



V,i\)+i\me^V,{\)+pS{-\) =0, 
^ Vl'\X) +tXm,^ Vl'^X) ± p Vl'\-X) = 0. 



Solutions of the Dirac equation are classified as follows 
1 

"V2 



x: 



A solution to Eq.(|26|) is obtained by introducing a function ax such that 

ax dfx 



S{X)=axfx, V,{-X) 



p dt 



dax 
~dt 



-i Xnie $ ax, 



and similarly for V^j'^'^(A), i.e. 



Here the Higgs background is rewritten as 



$ = iVfc cn{Hkt, k) = A^fc dn r. 



kM„t 



k 



M t, for k oo, 



using a relation that is based on the Jacobi's real transformation 

cn {Hkt, k) = dn (j, , 
dn {Hkt, k) = cn ^r, -^^ , 



sn {Hkt, k) = k sn ( r, — ) . 



Moreover, for k — > oo, the following approximations hold: 

cn ^r, — ^ ~ cos r + ~ 2 ^ ^' 

dn(r,i)~l-^sinV, 



sn ( T, — ^ ~ sin r — - (r sin 2r ) cos r. (32) 



With r = rrie/Mjj we obtain a solution for a\, 

axi^r) = exp |— 2zAr arcsin sn ^r, — ^ |, 



exp{-2 zArr + C j^i^ } = exp{-zAme t + C» (^^j }. (33) 



The function /, therefore, satisfies the following differential equation: 



where we have introduced a new parameter 



^'^^-4^Ardnfr,i') + g^A = 0, (34) 



if 

The function dn is doubly-periodic with periods 2 K and 4 i K'. Due to the periodicity we can 
discuss all properties of the elliptic function in the so-called fundamental period parallelogram 
which for dn is r = 2 ,^ K + 4 ir/ K' with < C,,i] < 1. An irreducible set of poles is given by 

/? = /3oo = ^K', /5' = /5oi=3zK', (36) 

with residues —i and +i respectively. Therefore, the singular points of the second order differential 
equation for fx are t = f3, f3' and their congruent points. The corresponding exponents are 0, 1 ± 
4 Ar for [3{[3'). We know that Eq. (|3^) possesses a fundamental set of solutions but, unfortunately, 
the exponents are not unequal integers and, therefore, we cannot apply the Hermite, Picard, 
Mittag-Leffler, Floquet theorem |p stating that the solutions are doubly-periodic functions of 
second kind and, in general, expressible as products of ratios of weierstrassian cr-functions, see 
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We have not been able to find an exact, explicit, solution to Eg . (|3^) but an approximated one 
can be given by using the WKBJ method, based on the observation that, for /c — >^ oo, 

<i> = -NkHkdn{Hkt,t) sn{Hkt,k) ~ cos r sin r, 

$ = iVfc cn {Hkt, ^) ~ 1 - ^ sin^ r. (37) 
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Therefore, in this hmit, we obtain 

/±(r)=exp{z duO^iu),}, 



e^(M) = 2Ardn(M, f ) ± g^ + 4rMn^(u 



1/2 



The integrals appearing in Eg . (|38|) give 

1 



du 



du dn M, 



arcsm sn r, 



k 



1 



+ 4 r^ dn^ 



1/2 



du 

2VnE I T. 



4F V" 2 

1/2 



- sin2r') + O (A; 



47^-/^2 sin^ M + O (A;-^) 
h 



2Vn 



2 V^r- - 
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r-i sin(2r)] 



where we have introduced the following parameters 

1 



-2,2 V^ + ml 



and where E' is the elliptic integral of second kind. In the limit A; — oo we also have 



n 



p 



M2 2 M2 fc2 ■ 



2VnT ~ St. 



To summarize, in WKBJ approximation, we find the following result for B of Eg. ([38|) , 



e^(u) = Ar 



2 — — sin^ u 



sin^ u^+O (k-^) . 



(38) 



(39) 



(40) 



(41) 



(42) 



Note that in the product ax fx (see Eg.(PB|)) arcsin(sn(r, 1/A;)) drops out and, by interacting with 
the time-dependent Higgs background, the positive(negative) energy plane wave-solutions that 
would correspond to a free electron of mass me. 



exp 



(43) 



receive a distortion factor which can be made arbitrarily small for large values of k. From Eg.(^) 
we derive 



S^{X) = axfx=exp{iJ^ duQ^iu),}, 



v,H-x) = ^sH^){x f (i 



sin^ r ' 
2P 



P 



mz 



A£^k^ 2£^k^^ 



sin^ r 



}■ 



(44) 



A similar results holds for the transverse components Vj|'^(A). A more accurate version of the 
approximated solution will now be derived. Starting from Eq.(0) we write 



/ = exp{9) F, 9 = 2 iXi am(r) 
with am(r) defined in Eq.(^ and where F is a solution of 



+ 



+ 4 r^ dn^ ( r, f ) - 2 sn ( r, f ) cn ( r, f 



1 



A;2 



1 



k 



0. 



(45) 



(46) 



The standard WKBJ solution of the above equation, based on the fact that $ is a slowly varying 
function of r, follows by introducing 



QKu) = + 4r^ - 4^ sin^M - 22A-^ sinMcosM + O (k~^) . 

1.2 h.2 V / 



We easily derive a solution for F, 



Qx{r) 

where the integral in the exponent gives 



exp|±i J duQx{u)^, 



(47) 



(4J 



1 r^ 

duQxiu) = 2VnT- — ^ 
y/n 



Note that for = we have an exact solution of Eq. 
S{\), Vl{X), or Vl{X). The corresponding equation becomes 

^ + 2ardn(r,i) /, = 0, 



T-- sin(2 t))+^^ (cos(2 t) - 1) + O (A;-^) . (49) 



Let fx be any of the functions 



sn T. 



with a solution 

fx = exp|— 2 iAr arcsm 

with am(r) defined in Eq.(^. Expanding for large values of k'^ gives 

, . 1 



— ^ I = exp|— 2 iAr am(r)|, 



am[T) = T 



4F 



r — sm T cos r 



(50) 
(51) 

(52) 



same conclusion by writing Eq.(|T^) as 

d^e 



Since 2 rr ~ rrie t we again recover the correct limit of the free Dirac equation. We arrive at the 

+ im-e $7^6 + mg$^ e = 0, (53) 
e= e'(t)w*, (54) 



dt' 



and expanding the spinor e{t) 



i=lA 
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where the u are eigenfunctions of the spin generator. It follows 

^ + 4r^d.(r.i)e'±2.^c.(r.i)s„(r.i)e = 0. (55) 

where the ± refers to i = 1, 2 and i = 3, 4 respectively. Using Eg. (^71) we easily derive, in WKBJ 
approximation, 

e^'^ = — ^— exp(-z /" duQ_{u)}, e^'^ = — ^— expj+i /" duQ+{u)}. (56) 
Q-{t) ^ Jo > Q+iV ^ JO 

In the above result it is understood that p = and \/TZ = r. 
4 Conclusions 

In this paper we have tried to answer an important question related to the possibility that the 
fundamental constants of nature might vary with time: can a time variation in the electron mass 
(any mass) be made formally consistent with the structure of the Standard Model of strong and 
electroweak interactions? For an exhaustive study of how a change in rrie alters physical effects 
we refer to 0. 

We have shown that the classical equation of motion for the Higgs field in the Minimal Standard 
Model admits a time- dependent solution which can be given in terms of the Jacobian elliptic 
function cn which is a doubly-periodic function of time. Therefore, we have a consistent picture 
of a time-dependent vev of the Higgs field. The background is not an external field but rather it 
is derived directly from the equations of motion. 

By choosing the primitive period of the Higgs vev large enough we have a picture of the vacuum 
that is not in any evident contradiction with plain experimental evidence. By this we mean that 
the present value of the proton-to-electron mass ratio is /i = 1836.1526645(57) and any significant 
variation of this parameter over a small time interval is excluded (A/i//i < 2 x 10"'' [0)- 

Next we have considered the coupling of fermions with the Higgs sector of the Minimal Standard 
Model and have studied the effect of their propagation in the time-dependent Higgs background. 

We have shown that the Dirac equation for the electron field is equivalent to a second order 
differential equation with doubly-periodic coefficients. In the limit of very large primitive period 
of the Higgs background this equation can be solved in WKBJ approximation. 

We summarize our findings by saying that the solutions to the Dirac equation in the slowly 
varying scalar-background field are plane-waves with a time-dependent distortion factor which 
becomes arbitrarily small in the limit of infinite primitive period of the Higgs vev. Therefore 
fermion fluctuations are almost like plane- waves with three- momentum p and mass mg but, strictly 
speaking, we have no time-dependent mass and £ of Eg. (^11) is not the energy of the state; we 
rather register a departure from the plane-wave shape through a time-dependent factor which, of 
course, could be fitted with a plane- wave solution with an effective me(t) parameter. 
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